Pairing symmetry of superconductivity coexisting with antiferromagnetism 
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Pairing symmetry in the superconducting state coexisting with antiferromagnetic order is studied 
based on a microscopic theory. We calculate the linearized Eliashberg's equation within the random 
phase approximation in the Hubbard model with a staggered field. We find that odd-frequency spin- 
triplet (equal-spin) s-wave pairing state can be realized. This result contradicts a naive expectation 
that antiferromagnetic order induces antiferromagnetic spin fiuctuation and favors spin-singlet d- 
wave pairing as in the standard strongly correlated systems. 



In condensed matter physics, coexistence of supercon- 
ductivity (SC) and magnetism has been a fundamental 
issue attracting a great deal of attention [ij. While SC 
coexists with ferromagnetism in, e.g., UGe2 and URhGe 
0,9; SC coexists with antiferromagnetism (AF) in, e.g., 
CeCu2Si2 and CeRhlns 0, Considering the fact that 
antiferromagnetic spin fluctuation mediates spin-singlet 
pairing, we simply expect that coexistence of antiferro- 
magnetic order (AFO) makes spin-singlet d-wave pairing 
more stable. However, pairing symmetry of Cooper pair 
in SC coexisting with AF is not simple. In CeCu2Si2 
and CeRhlns [6|-l8|, gapless excitation of quasiparticle 
has been reported experimentally. It has been pointed 
theoretically that a competition between even-frequency 
(even-w) spin-singlet d-wave and odd-frequency (odd-oj) 
[sl spin-singlet p-wave pairings is possible [l3l in the su- 
perconducting state due to the coexistence of AF. 

Odd-w pairing has been originally proposed by 
Berezinskii in 1974 A possibility of odd-o; spin-triplet 
pairing has been discussed in the context of superfluid 
•^He [9|. After that, Balatsky and Abrahams have pro- 
posed odd-o; spin-singlet pairing [ll| . Starting from these 
proposals, a lot of intensive studies about odd-o; pairing 
have been done both in bulk l^, 12-l3l and inhomoge- 
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systems. 

Although there have been several studies about SC co- 
existing with AFO, pairing symmetry of Cooper pair has 
not been fully resolved based on a microscopic theory. In 
the present paper, we study pairing symmetry in the su- 
perconducting state coexisting with AF focusing on the 
possible odd-oj spin-triplet pairing based on a microscopic 
theory. We consider the Hubbard model with a staggered 
field originating from commensurate AFO and solve the 
linearized Eliashberg's equation within the random phase 
approximation (RPA). We find that a competition be- 
tween even-oj spin-singlet d-wave and odd-o; spin-triplet 
(equal-spin) s-wave pairings is induced by the staggered 
field. We show that odd-w equal-spin-triplet s-wave pair- 
ing state becomes the most dominant under a high stag- 
gered field due to the charge fluctuation induced by AFO. 
This result contradicts a naive expectation that AFO in- 
duces antiferromagnetic spin fluctuation and favors spin- 
singlet d-wave pairing. 

We model a system where SC coexists with AF by sep- 



arately considering superconducting and antiferromag- 
netic electrons, i.e., we assume that interactions between 
superconducting and antiferromagnetic electrons are not 
so strong. We apply the Hubbard model to itinerant elec- 
trons generating superconductivity taking into account a 
background staggered field due to AFO by localized elec- 
trons as shown in Fig. [TJ We consider the nearest neigh- 
bor hopping t and the on-site Coulomb repulsion U and 
assume that the staggered field by AFO is a given quan- 
tity hs without solving dynamics of localized electrons. 
There are two kinds of sublattices A and B under finite 
hs. Potential felt by an itinerant electron with spin a is 
shifted by the staggered field +{—)schs on a sublattice A 
(B) with s-t-(_|,) = -|-(— )1, i.e., a sublattice A is J,-electron- 
rich while a sublattice B is ^-electron-rich under finite hs . 
The corresponding Hamiltonian is given by 
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FIG. 1: (Color online) The model considered in the present 
study with a hopping t, the on-site Coulomb repulsion U, and 
a staggered field ha. 
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where c'fj annihilates (creates) an itinerant electron with 
spin (T at a site i and number operator, — 1 

on a sublattice A while = — 1 on a sublattice B. The 
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Green's function without U is given by 
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where fc = (iw„,fc) is an elhpsis notation with the 
fermionic Matsubara frequency a;„ — {2n — 1)ttT and 
the momentum fe, /i is the chemical potential, e^^ is an 
(a, /3)-element of ffco-, / is a unit matrix, and a is the 
lattice constant of the unit cell. 

Within the RPA, we solve the linearized Eliashberg's 
equation 
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where A is an eigenvalue and N is the number of sites. 
A"^, (fc) and V^^,{k,k') are (a, /3)-elements of a gap 
function Ao-o-' (k) and an effective pairing interaction 
Vcrcr'{k,k'), respectively. T is equal to the supercon- 
ducting transition temperature Tq when A reaches unity. 
Thus, larger A corresponds to more stable superconduct- 
ing state. In Eq. dT)), a Cooper pair with equal (opposite) 
spins requires a' — cr(f7). V^a' [k, k') is given by 

k') =U U'xlfik k') U^tiik + fc'), (8) 
%.{k,k') = -U^xlf{k~k% (9) 

where U = UI and j(^{q) (x^^ (?)) denotes a part of 
susceptibility from bubble-type (ladder-type) diagrams. 
Here, q = (it'm, O) is an ellipsis notation with the bosonic 
Matsubara frequency Vm = ^mirT and the momentum q. 
xla^iq) and x]?g{q) are given by 
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where x^^^,{q) is the irreducible susceptibility. An (a,/3)- 
element of x^a' il) given by 
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X^^{q) and Xct^ (?) give the transverse spin, longitudinal 
spin, and charge susceptibilities written as 
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Note that the antiferromagnetic moment is assumed to 
be parallel to the spin quantization axis (z axis) in the 
present study. xtp~{q) and xfp('Z) indicate spin fluctua- 
tions in a xy plane and a z direction while Xch{q) indi- 
cates charge fluctuation. We define the maximum eigen- 
value of UxlsiqWXaaiq) (UXsaiq)) as the Stoner's fac- 
tor f^"^ for Xal^iq) ifs""^ for x^^ (<?)) and the larger one 
as the Stoner's factor fs- 

Hereafter, we consider Cooper pairs with following 
symmetries; even-w spin-singlet even-parity, even-oj spin- 
triplet odd-parity, odd-w spin-singlet odd-parity, and 
odd-w spin-triplet even-parity symmetries for various hg . 
Spin-singlet pairing has a spin state lit) (total pair 

spin Sz = 0) while spin-triplet pairing has spin states |tt) 
(5, = 1), lU) (S. = -1), and |n)+ lit) [S. = 0). We 
present that odd-oj spin-triplet {Sz = ±1) s-wave pair- 
ing state exceeds conventional even-cj spin-singlet d-wave 
one. In the actual calculation, fi in Eq. ([2]) is determined 
to tune 0.8 electrons per site. Temperature is chosen to 
be T = O.Oli. a;„ has a value from — (2nmax — 1)'!tT to 
(2n,„ax ~ l)7rr with rimax = 2048. Vrn has a value from 
— 2m,„ax7''T to 2TOniax7''T with TOmax = 2048. Wc take 
TV = 64 X 64. A gap function is normalized with a condi- 
tion J2k a 13 \ ^k)\'^ — 1' ^c choosc a valuc of U for /s 
to get 0.98 because /s determines a scale of {k,k'), 
to which A is sensitive. 

Fig. shows hs dependences of A. Here, we plot only 
pairings which become dominant. Odd-w spin-triplet 
(Sz = ±1) s-wave pairing state is the most stable one 
in a high hs region [hs > 2t) while even-w spin-singlet 
d-wave pairing state is the most stable one in a low hg 
region (0 < /is < 2t). Figs. [3] (a) and (b) show fc depen- 
dences of gap functions for even-w spin-singlet d-wave 
pairing at /is = and odd-w spin-triplet (5*2 — 1) s-wave 
one at hs = 2.6t, respectively, in the band basis. Here, 
we show only one of two bands where the Fermi surface 
exists for each pairing. The present odd-w spin-triplet 
(Sz = 1) s-wave pairing in a high hs region is on-site 
equal-spin one and has a weak modulation in real space. 

We discuss why the pairing competition between even- 
cj spin-singlet d-wave and odd-w spin-triplet (Sz = ±1) 
s-wave pairings occurs under hs as shown in Fig. [2j This 
pairing competition is caused by suppression of spin fluc- 
tuation in a xy plane and enhancement of charge fluc- 
tuation originating from hs. xt^iq), Xl^{q), and Xch(g) 
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FIG. 2: (Color online) hs dependences of A for leading pairing 
symmetries; even-oj spin-singlet d-wave and odd-oj spin-triplet 
(5z — ±1) s-wave. 




FIG. 4: (Color online) hs dependences of xt^{q), xlpiq), 
and Xch('?). The maximum values after diagonalization are 
plotted. 
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FIG. 3: (Color online) fc dependences of gap functions for (a) 
even-oj spin-singlet d-wave pairing aX hs = Q and (b) odd-o; 
spin-triplet {Sz = 1) s-wave one at ha = 2.Gt in the band 
basis at aj„ = ttT. Dotted (green) and solid (black) lines 
denote nodal lines and the Fermi surfaces, respectively. 
and "— " denote signs of gap functions. 



depend on hg as shown in Fig. 21 With increasing hg, 
Xsp~(9) decreases while xl^{q) and Xch(9) increase. In 
other words, spin fluctuation in a xy plane is suppressed 
while that in a ^ direction and charge fluctuation are en- 
hanced. In order to understand how these changes influ- 
ence the symmetry of Cooper pair, we look at the single- 
band Hubbard model without hg . In the case where spin 
or charge fluctuation dominates, an effective pairing in- 
teraction Vaa'{k,k') in the single-band Hubbard model 
without hs is written as 

~ ]-U^ [2x+p-(fc + fc') + Xi;{k - k') - xch(fc - fc')] , 



(17) 
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with the transverse spin susceptibility (<?) , the longi- 
tudinal spin susceptibility xlp il) i and the charge suscep- 
tibility Xch('?)j which are positive numbers. As is shown 
in Eqs. ([T7)) and suppression of the spin fluctuation 
in a xy plane destabilizes pairing with Sz = while en- 
hancement of charge fluctuation destabilizes (stabilizes) 
pairing with = {Sz — ±1). By these features, 
hs causes the pairing competition between even-w spin- 
singlet d-wave and odd-w spin-triplet (Sz = ±1) s-wave 
pairings. 

Here, you may wonder why not even-w spin-triplet 
{Sz = ±1) p-wave but odd-w spin-triplet {Sz = ±1) s- 
wave pairing is stabilized. The reason is that an electron 
with each spin feels inhomogeneity in real space due to 
the background AFO. Even-w spin-triplet {Sz — ±1) p- 
wave pairing tends to be equal-spin pairing between the 
nearest sites. However, in a high hs region, t-electron- 
rich and 4,-electron-rich sites are adjacent. Therefore, 
x!^"^'""(9), which mediates equal-spin pairing between 
the nearest sites, is suppressed while x\i^'^^{l) ^^^d 
x'j!"^'^^((7), which mediate on-site equal-spin pairing, 
are enhanced (see Fig. [5]). Thus even-w spin-triplet 
{Sz = ±1) p-wave pairing is not allowed there and odd-o; 
spin-triplet {Sz = ±1) s-wave one appears. 

Next, we discuss why xi~(g) decreases while Xsp('?) 
and Xch{<l) increase with increasing hs as shown in Fig. 
m The decrease of xtp'il) ™d the increase of Xspil) 
can be interpreted as the suppression of spin fluctua- 
tion in a xy plane and the enhancement of that in a 
z direction, respectively, due to hs along a z direction. 
The increase of Xchl?) is caused by inhomogeneity in real 
space for an electron with each spin due to hs- x^^{q), 
which gives Xspil) ^^'^ Xch{<l) as written in Eqs. ([T5|) and 
(fT6l) . depends on hs as shown in Fig. [5l Xct"^'"^('7) and 
X^^s''"^{q) ars opposite in sign to each other. In a low hs 
region, they have comparable absolute values and cancel 
out each other in the summation of right side of Xcif (?) 
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FIG. 5: (Color online) dependences of (AA) xl'^J^'^'^{<l), 

(AB) x^'^^C?), (BA) x^".'?'^*(a), and (BB) xl'^f'^'^iq)- In 
each panel, the values at g-point where the absolute values 
get maximum are plotted for all four combinations of spins. 



Eq. (IT51) , while they enhance each other in that of right 
side of Xsp'"^(.q) Eq. (|T5|) . This picture is similar to 
that in the single-band Hubbard model without hs- In 
a high hs region, on the other hand, the above picture 
breaks down since only the susceptibility between a spe- 
cific combination of spins is prominently enhanced. As 
shown in Fig. [5l the magnitudes of xli^'^^il) in ^ panel 
(AA), x\^^'^^{q) in a panel (AB), x\i^'^^{<l) in a panel 
(BA), and X^^^'^^il) in a panel (BB) are prominently 
enhanced. This feature arises from the inhomogeneity in 
real space for an electron with each spin, i.e., 4,-electron- 
richness on a sublattice A and f-electron-richness on a 
sublattice B due to the presence of hg. In the present 
situation, the above cancellation for low hs in the ex- 
pression of Xch{q) does not work any more. Thus, in a 
high hs region, Xchil) becomes dominant in addition to 

In summary, we have studied pairing symmetry in the 
superconducting state coexisting with AF by solving the 
linearized Eliashberg's equation in the Hubbard model 
with a staggered field using the RPA. As a result, we 
have found that odd-w equal-spin-triplet s-wave pairing 
can be realized. This is caused by the suppression of spin 
fiuctuation in a plane perpendicular to an antiferromag- 
nctic moment and the enhancement of charge fiuctuation 
due to the background AFO. In particular, the enhance- 



ment of charge fluctuation originates from inhomogeneity 
in real space for an electron with each spin induced by the 
background AFO. The present result contradicts a naive 
picture that AFO might help antiferromagnetic spin fluc- 
tuation and favor even-w spin-singlet d-wave pairing. We 
hope that the odd-w spin-triplet pairing will be verified 
experimentally in strongly correlated systems where SC 
coexists with AF. 
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